Structural insights reveal interplay between LAG-3 homodimerization, ligand binding, and function.

Lymphocyte activation gene-3 (LAG-3) is an inhibitory receptor expressed on activated T cells and an emerging immunotherapy target. Domain 1 (D1) of LAG-3, which has been purported to directly interact with major histocompatibility complex class II (MHCII) and fibrinogen-like protein 1 (FGL1), has been the major focus for the development of therapeutic antibodies that inhibit LAG-3 receptor-ligand interactions and restore T cell function. Here, we present a high-resolution structure of glycosylated mouse LAG-3 ectodomain, identifying that cis-homodimerization, mediated through a network of hydrophobic residues within domain 2 (D2), is critically required for LAG-3 function. Additionally, we found a previously unidentified key protein-glycan interaction in the dimer interface that affects the spatial orientation of the neighboring D1 domain. Mutation of LAG-3 D2 residues reduced dimer formation, dramatically abolished LAG-3 binding to both MHCII and FGL1 ligands, and consequentially inhibited the role of LAG-3 in suppressing T cell responses. Intriguingly, we showed that antibodies directed against D1, D2, and D3 domains are all capable of blocking LAG-3 dimer formation and MHCII and FGL-1 ligand binding, suggesting a potential allosteric model of LAG-3 function tightly regulated by dimerization. Furthermore, our work reveals unique epitopes, in addition to D1, that can be targeted for immunotherapy of cancer and other human diseases.

Deep residual networks for crystallography trained on synthetic data.

The use of artificial intelligence to process diffraction images is challenged by the need to assemble large and precisely designed training data sets. To address this, a codebase called Resonet was developed for synthesizing diffraction data and training residual neural networks on these data. Here, two per-pattern capabilities of Resonet are demonstrated: (i) interpretation of crystal resolution and (ii) identification of overlapping lattices. Resonet was tested across a compilation of diffraction images from synchrotron experiments and X-ray free-electron laser experiments. Crucially, these models readily execute on graphics processing units and can thus significantly outperform conventional algorithms. While Resonet is currently utilized to provide real-time feedback for macromolecular crystallography users at the Stanford Synchrotron Radiation Lightsource, its simple Python-based interface makes it easy to embed in other processing frameworks. This work highlights the utility of physics-based simulation for training deep neural networks and lays the groundwork for the development of additional models to enhance diffraction collection and analysis.

Remote laboratory training for high school students: grocery store based hands-on project in protein crystallography.

The COVID-19 pandemic measures forced students to stay home and confined them to remote learning. This had a large impact on laboratory experiments, which are often impossible to complete from home. This article is a resource for instructors/educators to introduce the topic of structural biology and crystallographic methods. The main focus is to describe a hands-on crystallization laboratory exercise that can be carried out remotely at home with safe household products. X-ray crystallography is a vital technique for determining protein structure and function. This information can be used to understand fundamental biological processes and to help in the design of life-saving medications. Here, a method was developed to teach crystallography using reagents and equipment that can be found in grocery stores. The steps involved in a crystallography experiment are detailed with links and references to additional resources.

Triepoxide formation by a flavin-dependent monooxygenase in monensin biosynthesis.

Monensin A is a prototypical natural polyether polyketide antibiotic. It acts by binding a metal cation and facilitating its transport across the cell membrane. Biosynthesis of monensin A involves construction of a polyene polyketide backbone, subsequent epoxidation of the alkenes, and, lastly, formation of cyclic ethers via epoxide-opening cyclization. MonCI, a flavin-dependent monooxygenase, is thought to transform all three alkenes in the intermediate polyketide premonensin A into epoxides. Our crystallographic study has revealed that MonCI's exquisite stereocontrol is due to the preorganization of the active site residues which allows only one specific face of the alkene to approach the reactive C(4a)-hydroperoxyflavin moiety. Furthermore, MonCI has an unusually large substrate-binding cavity that can accommodate premonensin A in an extended or folded conformation which allows any of the three alkenes to be placed next to C(4a)-hydroperoxyflavin. MonCI, with its ability to perform multiple epoxidations on the same substrate in a stereospecific manner, demonstrates the extraordinary versatility of the flavin-dependent monooxygenase family of enzymes.

Enhanced active-site electric field accelerates enzyme catalysis.

The design and improvement of enzymes based on physical principles remain challenging. Here we demonstrate that the principle of electrostatic catalysis can be leveraged to substantially improve a natural enzyme's activity. We enhanced the active-site electric field in horse liver alcohol dehydrogenase by replacing the serine hydrogen-bond donor with threonine and replacing the catalytic Zn2+ with Co2+. Based on the electric field enhancement, we make a quantitative prediction of rate acceleration-50-fold faster than the wild-type enzyme-which was in close agreement with experimental measurements. The effects of the hydrogen bonding and metal coordination, two distinct chemical forces, are described by a unified physical quantity-electric field, which is quantitative, and shown here to be additive and predictive. These results suggest a new design paradigm for both biological and non-biological catalysts.

Discovery and Characterization of Antibody Probes of Module 2 of the 6-Deoxyerythronolide B Synthase.

Fragment antigen-binding domains of antibodies (Fabs) are powerful probes of structure-function relationships of assembly line polyketide synthases (PKSs). We report the discovery and characterization of Fabs interrogating the structure and function of the ketosynthase-acyltransferase (KS-AT) core of Module 2 of the 6-deoxyerythronolide B synthase (DEBS). Two Fabs (AC2 and BB1) were identified to potently inhibit the catalytic activity of Module 2. Both AC2 and BB1 were found to modulate ACP-mediated reactions catalyzed by this module, albeit by distinct mechanisms. AC2 primarily affects the rate (kcat), whereas BB1 increases the KM of an ACP-mediated reaction. A third Fab, AA5, binds to the KS-AT fragment of DEBS Module 2 without altering either parameter; it is phenotypically reminiscent of a previously characterized Fab, 1B2, shown to principally recognize the N-terminal helical docking domain of DEBS Module 3. Crystal structures of AA5 and 1B2 bound to the KS-AT fragment of Module 2 were solved to 2.70 and 2.65 Å resolution, respectively, and revealed entirely distinct recognition features of the two antibodies. The new tools and insights reported here pave the way toward advancing our understanding of the structure-function relationships of DEBS Module 2, arguably the most well-studied module of an assembly line PKS.

Structural and functional analyses of the echinomycin resistance conferring protein Ecm16 from Streptomyces lasalocidi.

Echinomycin is a natural product DNA bisintercalator antibiotic. The echinomycin biosynthetic gene cluster in Streptomyces lasalocidi includes a gene encoding the self-resistance protein Ecm16. Here, we present the 2.0 Å resolution crystal structure of Ecm16 bound to adenosine diphosphate. The structure of Ecm16 closely resembles that of UvrA, the DNA damage sensor component of the prokaryotic nucleotide excision repair system, but Ecm16 lacks the UvrB-binding domain and its associated zinc-binding module found in UvrA. Mutagenesis study revealed that the insertion domain of Ecm16 is required for DNA binding. Furthermore, the specific amino acid sequence of the insertion domain allows Ecm16 to distinguish echinomycin-bound DNA from normal DNA and link substrate binding to ATP hydrolysis activity. Expression of ecm16 in the heterologous host Brevibacillus choshinensis conferred resistance against echinomycin and other quinomycin antibiotics, including thiocoraline, quinaldopeptin, and sandramycin. Our study provides new insight into how the producers of DNA bisintercalator antibiotics fend off the toxic compounds that they produce.

Potent and selective covalent inhibition of the papain-like protease from SARS-CoV-2.

Direct-acting antivirals are needed to combat coronavirus disease 2019 (COVID-19), which is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The papain-like protease (PLpro) domain of Nsp3 from SARS-CoV-2 is essential for viral replication. In addition, PLpro dysregulates the host immune response by cleaving ubiquitin and interferon-stimulated gene 15 protein from host proteins. As a result, PLpro is a promising target for inhibition by small-molecule therapeutics. Here we design a series of covalent inhibitors by introducing a peptidomimetic linker and reactive electrophile onto analogs of the noncovalent PLpro inhibitor GRL0617. The most potent compound inhibits PLpro with kinact/KI = 9,600 M-1 s-1, achieves sub-µM EC50 values against three SARS-CoV-2 variants in mammalian cell lines, and does not inhibit a panel of human deubiquitinases (DUBs) at >30 µM concentrations of inhibitor. An X-ray co-crystal structure of the compound bound to PLpro validates our design strategy and establishes the molecular basis for covalent inhibition and selectivity against structurally similar human DUBs. These findings present an opportunity for further development of covalent PLpro inhibitors.

Protein Electric Fields Enable Faster and Longer-Lasting Covalent Inhibition of ß-Lactamases.

The widespread design of covalent drugs has focused on crafting reactive groups of proper electrophilicity and positioning toward targeted amino-acid nucleophiles. We found that environmental electric fields projected onto a reactive chemical bond, an overlooked design element, play essential roles in the covalent inhibition of TEM-1 ß-lactamase by avibactam. Using the vibrational Stark effect, the magnitudes of the electric fields that are exerted by TEM active sites onto avibactam's reactive C═O were measured and demonstrate an electrostatic gating effect that promotes bond formation yet relatively suppresses the reverse dissociation. These results suggest new principles of covalent drug design and off-target site prediction. Unlike shape and electrostatic complementary which address binding constants, electrostatic catalysis drives reaction rates, essential for covalent inhibition, and deepens our understanding of chemical reactivity, selectivity, and stability in complex systems.

Structure-Based Prototyping of Allosteric Inhibitors of Human Uridine/Cytidine Kinase 2 (UCK2).

Pyrimidine nucleotide biosynthesis in humans is a promising chemotherapeutic target for infectious diseases caused by RNA viruses. Because mammalian cells derive pyrimidine ribonucleotides through a combination of de novo biosynthesis and salvage, combined inhibition of dihydroorotate dehydrogenase (DHODH; the first committed step in de novo pyrimidine nucleotide biosynthesis) and uridine/cytidine kinase 2 (UCK2; the first step in salvage of exogenous nucleosides) strongly attenuates viral replication in infected cells. However, while several pharmacologically promising inhibitors of human DHODH are known, to date there are no reports of medicinally viable leads against UCK2. Here, we use structure-based drug prototyping to identify two classes of promising leads that noncompetitively inhibit UCK2 activity. In the process, we have identified a hitherto unknown allosteric site at the intersubunit interface of this homotetrameric enzyme. By reducing the kcat of human UCK2 without altering its KM, these new inhibitors have the potential to enable systematic dialing of the fractional inhibition of pyrimidine salvage to achieve the desired antiviral effect with minimal host toxicity.

Structural and functional characterization of NEMO cleavage by SARS-CoV-2 3CLpro.

In addition to its essential role in viral polyprotein processing, the SARS-CoV-2 3C-like protease (3CLpro) can cleave human immune signaling proteins, like NF-κB Essential Modulator (NEMO) and deregulate the host immune response. Here, in vitro assays show that SARS-CoV-2 3CLpro cleaves NEMO with fine-tuned efficiency. Analysis of the 2.50 Å resolution crystal structure of 3CLpro C145S bound to NEMO226-234 reveals subsites that tolerate a range of viral and host substrates through main chain hydrogen bonds while also enforcing specificity using side chain hydrogen bonds and hydrophobic contacts. Machine learning- and physics-based computational methods predict that variation in key binding residues of 3CLpro-NEMO helps explain the high fitness of SARS-CoV-2 in humans. We posit that cleavage of NEMO is an important piece of information to be accounted for, in the pathology of COVID-19.

Protein sequence design with a learned potential.

The task of protein sequence design is central to nearly all rational protein engineering problems, and enormous effort has gone into the development of energy functions to guide design. Here, we investigate the capability of a deep neural network model to automate design of sequences onto protein backbones, having learned directly from crystal structure data and without any human-specified priors. The model generalizes to native topologies not seen during training, producing experimentally stable designs. We evaluate the generalizability of our method to a de novo TIM-barrel scaffold. The model produces novel sequences, and high-resolution crystal structures of two designs show excellent agreement with in silico models. Our findings demonstrate the tractability of an entirely learned method for protein sequence design.

Energetic Basis and Design of Enzyme Function Demonstrated Using GFP, an Excited-State Enzyme.

The past decades have witnessed an explosion of de novo protein designs with a remarkable range of scaffolds. It remains challenging, however, to design catalytic functions that are competitive with naturally occurring counterparts as well as biomimetic or nonbiological catalysts. Although directed evolution often offers efficient solutions, the fitness landscape remains opaque. Green fluorescent protein (GFP), which has revolutionized biological imaging and assays, is one of the most redesigned proteins. While not an enzyme in the conventional sense, GFPs feature competing excited-state decay pathways with the same steric and electrostatic origins as conventional ground-state catalysts, and they exert exquisite control over multiple reaction outcomes through the same principles. Thus, GFP is an "excited-state enzyme". Herein we show that rationally designed mutants and hybrids that contain environmental mutations and substituted chromophores provide the basis for a quantitative model and prediction that describes the influence of sterics and electrostatics on excited-state catalysis of GFPs. As both perturbations can selectively bias photoisomerization pathways, GFPs with fluorescence quantum yields (FQYs) and photoswitching characteristics tailored for specific applications could be predicted and then demonstrated. The underlying energetic landscape, readily accessible via spectroscopy for GFPs, offers an important missing link in the design of protein function that is generalizable to catalyst design.

Stabilization of glucose-6-phosphate dehydrogenase oligomers enhances catalytic activity and stability of clinical variants.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a genetic trait that can cause hemolytic anemia. To date, over 150 nonsynonymous mutations have been identified in G6PD, with pathogenic mutations clustering near the dimer and/or tetramer interface and the allosteric NADP+-binding site. Recently, our lab identified a small molecule that activates G6PD variants by stabilizing the allosteric NADP+ and dimer complex, suggesting therapeutics that target these regions may improve structural defects. Here, we elucidated the connection between allosteric NADP+ binding, oligomerization, and pathogenicity to determine whether oligomer stabilization can be used as a therapeutic strategy for G6PD deficiency (G6PDdef). We first solved the crystal structure for G6PDK403Q, a mutant that mimics the physiological acetylation of wild-type G6PD in erythrocytes and demonstrated that loss of allosteric NADP+ binding induces conformational changes in the dimer. These structural changes prevent tetramerization, are unique to Class I variants (the most severe form of G6PDdef), and cause the deactivation and destabilization of G6PD. We also introduced nonnative cysteines at the oligomer interfaces and found that the tetramer complex is more catalytically active and stable than the dimer. Furthermore, stabilizing the dimer and tetramer improved protein stability in clinical variants, regardless of clinical classification, with tetramerization also improving the activity of G6PDK403Q and Class I variants. These findings were validated using enzyme activity and thermostability assays, analytical size-exclusion chromatography (SEC), and SEC coupled with small-angle X-ray scattering (SEC-SAXS). Taken together, our findings suggest a potential therapeutic strategy for G6PDdef and provide a foundation for future drug discovery efforts.

Photosynthetic reaction center variants made via genetic code expansion show Tyr at M210 tunes the initial electron transfer mechanism.

Photosynthetic reaction centers (RCs) from Rhodobacter sphaeroides were engineered to vary the electronic properties of a key tyrosine (M210) close to an essential electron transfer component via its replacement with site-specific, genetically encoded noncanonical amino acid tyrosine analogs. High fidelity of noncanonical amino acid incorporation was verified with mass spectrometry and X-ray crystallography and demonstrated that RC variants exhibit no significant structural alterations relative to wild type (WT). Ultrafast transient absorption spectroscopy indicates the excited primary electron donor, P*, decays via a ∼4-ps and a ∼20-ps population to produce the charge-separated state P+HA- in all variants. Global analysis indicates that in the ∼4-ps population, P+HA- forms through a two-step process, P*→ P+BA-→ P+HA-, while in the ∼20-ps population, it forms via a one-step P* → P+HA- superexchange mechanism. The percentage of the P* population that decays via the superexchange route varies from ∼25 to ∼45% among variants, while in WT, this percentage is ∼15%. Increases in the P* population that decays via superexchange correlate with increases in the free energy of the P+BA- intermediate caused by a given M210 tyrosine analog. This was experimentally estimated through resonance Stark spectroscopy, redox titrations, and near-infrared absorption measurements. As the most energetically perturbative variant, 3-nitrotyrosine at M210 creates an ∼110-meV increase in the free energy of P+BA- along with a dramatic diminution of the 1,030-nm transient absorption band indicative of P+BA- formation. Collectively, this work indicates the tyrosine at M210 tunes the mechanism of primary electron transfer in the RC.

Modular polyketide synthase contains two reaction chambers that operate asynchronously.

Type I modular polyketide synthases are homodimeric multidomain assembly line enzymes that synthesize a variety of polyketide natural products by performing polyketide chain extension and ß-keto group modification reactions. We determined the 2.4-angstrom-resolution x-ray crystal structure and the 3.1-angstrom-resolution cryo­electron microscopy structure of the Lsd14 polyketide synthase, stalled at the transacylation and condensation steps, respectively. These structures revealed how the constituent domains are positioned relative to each other, how they rearrange depending on the step in the reaction cycle, and the specific interactions formed between the domains. Like the evolutionarily related mammalian fatty acid synthase, Lsd14 contains two reaction chambers, but only one chamber in Lsd14 has the full complement of catalytic domains, indicating that only one chamber produces the polyketide product at any given time.

Reproducibility of protein x-ray diffuse scattering and potential utility for modeling atomic displacement parameters.

Protein structure and dynamics can be probed using x-ray crystallography. Whereas the Bragg peaks are only sensitive to the average unit-cell electron density, the signal between the Bragg peaks-diffuse scattering-is sensitive to spatial correlations in electron-density variations. Although diffuse scattering contains valuable information about protein dynamics, the diffuse signal is more difficult to isolate from the background compared to the Bragg signal, and the reproducibility of diffuse signal is not yet well understood. We present a systematic study of the reproducibility of diffuse scattering from isocyanide hydratase in three different protein forms. Both replicate diffuse datasets and datasets obtained from different mutants were similar in pairwise comparisons (Pearson correlation coefficient ≥0.8). The data were processed in a manner inspired by previously published methods using custom software with modular design, enabling us to perform an analysis of various data processing choices to determine how to obtain the highest quality data as assessed using unbiased measures of symmetry and reproducibility. The diffuse data were then used to characterize atomic mobility using a liquid-like motions (LLM) model. This characterization was able to discriminate between distinct anisotropic atomic displacement parameter (ADP) models arising from different anisotropic scaling choices that agreed comparably with the Bragg data. Our results emphasize the importance of data reproducibility as a model-free measure of diffuse data quality, illustrate the ability of LLM analysis of diffuse scattering to select among alternative ADP models, and offer insights into the design of successful diffuse scattering experiments.

Antibody Probes of Module 1 of the 6-Deoxyerythronolide B Synthase Reveal an Extended Conformation During Ketoreduction.

The 6-deoxyerythronolide B synthase (DEBS) is a prototypical assembly line polyketide synthase (PKS) that synthesizes the macrocyclic core of the antibiotic erythromycin. Each of its six multidomain modules presumably sample distinct conformations, as biosynthetic intermediates tethered to their acyl carrier proteins interact with multiple active sites during the courses of their catalytic cycles. The spatiotemporal details underlying these protein dynamics remain elusive. Here, we investigate one aspect of this conformational flexibility using two domain-specific monoclonal antibody fragments (Fabs) isolated from a very large naïve human antibody library. Both Fabs, designated 1D10 and 2G10, were bound specifically and with high affinity to the ketoreductase domain of DEBS module 1 (KR1). Comparative kinetic analysis of stand-alone KR1 as well as a truncated bimodular derivative of DEBS revealed that 1D10 inhibited KR1 activity whereas 2G10 did not. Co-crystal structures of each KR1-Fab complex provided a mechanistic rationale for this difference. A hybrid PKS module harboring KR1 was engineered, whose individual catalytic domains have been crystallographically characterized at high resolution. Size exclusion chromatography coupled to small-angle X-ray scattering (SEC-SAXS) of this hybrid module bound to 1D10 provided further support for the catalytic relevance of the "extended" model of a PKS module. Our findings reinforce the power of monoclonal antibodies as tools to interrogate structure-function relationships of assembly line PKSs.

A CLC-ec1 mutant reveals global conformational change and suggests a unifying mechanism for the CLC Cl-/H+ transport cycle.

Among coupled exchangers, CLCs uniquely catalyze the exchange of oppositely charged ions (Cl- for H+). Transport-cycle models to describe and explain this unusual mechanism have been proposed based on known CLC structures. While the proposed models harmonize with many experimental findings, gaps and inconsistencies in our understanding have remained. One limitation has been that global conformational change - which occurs in all conventional transporter mechanisms - has not been observed in any high-resolution structure. Here, we describe the 2.6 Å structure of a CLC mutant designed to mimic the fully H+-loaded transporter. This structure reveals a global conformational change to improve accessibility for the Cl- substrate from the extracellular side and new conformations for two key glutamate residues. Together with DEER measurements, MD simulations, and functional studies, this new structure provides evidence for a unified model of H+/Cl- transport that reconciles existing data on all CLC-type proteins.

Cells are shielded from harmful molecules and other threats by a thin, flexible layer called the membrane. However, this barrier also prevents chloride, sodium, protons and other ions from moving in or out of the cell. Channels and transporters are two types of membrane proteins that form passageways for these charged particles. Channels let ions flow freely from one side of the membrane to the other. To do so, these proteins change their three-dimensional shape to open or close as needed. On the other hand, transporters actively pump ions across the membrane to allow the charged particles to accumulate on one side. The shape changes needed for that type of movement are different: the transporters have to open a passageway on one side of the membrane while closing it on the other side, alternating openings to one side or the other. In general, channels and transporters are not related to each other, but one exception is a group called CLCs proteins. Present in many organisms, this family contains a mixture of channels and transporters. For example, humans have nine CLC proteins: four are channels that allow chloride ions in and out, and five are 'exchange transporters' that make protons and chloride ions cross the membrane in opposite directions. These proteins let one type of charged particle move freely across the membrane, which generates energy that the transporter then uses to actively pump the other ion in the direction needed by the cell. Yet, the exact three-dimensional changes required for CLC transporters and channels to perform their roles are still unknown. To investigate this question, Chavan, Cheng et al. harnessed a technique called X-ray crystallography, which allows scientists to look at biological molecules at the level of the atom. This was paired with other methods to examine a CLC mutant that adopts the shape of a normal CLC transporter when it is loaded with a proton. The experiments revealed how various elements in the transporter move relative to each other to adopt a structure that allows protons and chloride ions to enter the protein from opposite sides of the membrane, using separate pathways. While obtained on a bacterial CLC, these results can be applied to other CLC channels and transporters (including those in humans), shedding light on how this family transports charged particles across membranes. From bone diseases to certain types of seizures, many human conditions are associated with poorly functioning CLCs. Understanding the way these structures change their shapes to perform their roles could help to design new therapies for these health problems.

Electrostatic control of photoisomerization pathways in proteins.

Rotation around a specific bond after photoexcitation is central to vision and emerging opportunities in optogenetics, super-resolution microscopy, and photoactive molecular devices. Competing roles for steric and electrostatic effects that govern bond-specific photoisomerization have been widely discussed, the latter originating from chromophore charge transfer upon excitation. We systematically altered the electrostatic properties of the green fluorescent protein chromophore in a photoswitchable variant, Dronpa2, using amber suppression to introduce electron-donating and electron-withdrawing groups to the phenolate ring. Through analysis of the absorption (color), fluorescence quantum yield, and energy barriers to ground- and excited-state isomerization, we evaluate the contributions of sterics and electrostatics quantitatively and demonstrate how electrostatic effects bias the pathway of chromophore photoisomerization, leading to a generalized framework to guide protein design.